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Abstract 

We explore the possibility that the diffuse gamma-ray background radiation (GBR) 
at high galactic latitudes could be dominated by inverse Compton scattering of 
cosmic ray (CR) electrons on the cosmic microwave background radiation and on 
starlight from our own galaxy. Assuming that the mechanisms accelerating galactic 
CR hadrons and electrons are the same, we derive simple and successful relations 
between the spectral indices of the GBR above a few MeV, and of the CR electrons 
and CR nuclei above a few GeV. We reproduce the observed intensity and angular 
dependence of the GBR, in directions away from the galactic disk and centre, without 
recourse to hypothetical extragalactic sources. 



PACS numbers: 98.70.Sa, 98.70.Rz, 98.70.Vc. 



1 Introduction 



The existence of an isotropic, diffuse gamma background radiation (GBR) was first 
suggested by data from the SAS 2 satellite (Thompson & Fichtel 1982). The EGRET 
instrument on the Compton Gamma Ray Observatory confirmed this finding: by 
removal of point sources and of the galactic-disk and galactic-centre emission, and 
after an extrapolation to zero local column density, a uniformly distributed GBR 
was found, of alleged extragalactic origin (Sreekumar et al. 1998). Above an energy 
of ~ 10 MeV, this radiation -to which we shall refer throughout simply as "the 
GBR"- has a featureless spectrum, shown in Fig. 1, which is very well described 
by a simple power-law form, dF/dE oc E _/3 , with (3 2.10 ± 0.03 (Sreekumar et al. 
1998). 

The origin of the GBR is still unknown. The published candidate sources range 
from the quite conventional to the decisively speculative. Perhaps the most conser- 
vative hypothesis for the origin of an isotropic GBR is that it is extragalactic, and 
originates from active galaxies (Bignami et al. 1979; Kazanas & Protheroe 1983; 
Stecker & Salamon 1996). The fact that blazars have a 7-ray spectrum with an av- 
erage index 2.15 ± 0.04, compatible with that of the GBR, supports this hypothesis 
(Chiang & Mukerjee 1998). The possibility has also been discussed that Geminga- 
type pulsars, expelled into the galactic halo by asymmetric supernova explosions, be 
abundant enough to explain the GBR (Dixon et al. 1998; Hartmann 1995). More 
exotic hypotheses include a baryon-symmetric universe (Stecker et al. 1971), now 
excluded (Cohen et al. 1998), primordial black hole evaporation (Page & Hawking 
1976; Hawking 1977), supermassive black holes formed at very high redshift (Gnedin 
& Ostriker 1992), annihilation of weakly interactive big-bang remnants (Silk & Sred- 
nicki 1984; Rudaz & Stecker 1991), and a long etc. 

However, the EGRET GBR data in directions above the galactic disk and cen- 
tre show a significant deviation from isotropy, correlated with the structure of our 
galaxy and our position relative to its centre (Dar et al. 1999). This advocates a 
local (as opposed to cosmological) origin for the GBR. Indications of a large galac- 
tic contribution to the GBR at large latitudes were independently found by Dixon 
et al. (1998) by means of a wavelet-based "non-parametric" approach that makes 
no reference to a particular model. Strong & Moskalenko (1998) and Moskalenko 
& Strong (2000) also found that the contribution of inverse Compton scattering of 
galactic cosmic ray electrons to the diffuse 7-ray background is presumably much 
larger than previously thought. In this paper we go one step further and explore in 
detail the possibility (Dar et al. 1999) that the diffuse gamma-ray background radi- 
ation at high galactic latitudes could be dominated by inverse Compton scattering 
of cosmic ray (CR) electrons on the cosmic microwave background radiation and on 
starlight from our own galaxy. In Section 2 we briefly review the GBR data and the 
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evidence for its correlation with our position in the Galaxy. 

The CR-proton and CR-electron spectra are briefly reviewed in Section 3. The 
origin, spectrum and composition of non-solar cosmic ray protons and nuclei have 
been debated for almost a century The measurements now extend over some 30 
orders of magnitude in flux and some 15 orders of magnitude in energy, up to an 
astonishing E ~ 3 x 10 11 GeV (Bird et al. 1995, Takeda et al. 1998, Berezinskii et 
al. 1990 and references therein). Above ~5 GeV, this spectrum has also a power-law 
form E _/3 , with two small variations in the "index" (3 at the so-called "CR knee" and 
"CR ankle". The local spectrum of CR electrons, shown in Fig. 2, is much harder 
to measure; it is only known up to ~ 10 3 GeV and, above ~ 5 GeV, it is also well 
described by a simple power law. 

In Sections 4 and 5 we discuss relations between the indices of the GBR and 
the CR electron and proton spectra. In so doing, we make few and very simple 
assumptions: that the mechanism accelerating CR hadrons and CR electrons is the 
same (a moving magnetic "mirror"), that the locally- measured electron spectrum is 
representative of its average form throughout the Galaxy, that above a certain en- 
ergy, inevitably, the electron spectrum is modulated by inverse Compton scattering 
on starlight and on the microwave background radiation, and that the GBR is dom- 
inated by the resulting Compton up-scattered photons. This allows one to derive, 
successfully, the GBR index from the electron index and the electron index from the 
proton index. The GBR index, as observed by EGRET, is uncannily directionally 
uniform. We interpret this fact as strong support for our simple assumptions. 

In Section 6 we tackle a more difficult and potentially controversial subject: the 
origin and magnitude of the GBR. In a sense, our proposed explanation -that the 
GBR originates from inverse Compton scattering in our own galaxy (Dar et al. 1999) 
- is more conservative than any of the previously suggested origins. 

The non-conventional aspect of our hypothesis is that, in order to reproduce the 
observed intensity of the GBR, we must assume the scale height of our galaxy's 
CR-electron distribution to be almost twice the traditionally-accepted upper limit. 
Because of this, in Section 6, we briefly review the basis of the conventional wisdom 
and our critical view of it, whose main points are the following. Moskalenko, Strong 
and their collaborators have developed a very detailed understanding of the CR, 
radio and 7 observations of our galaxy. To fit the data, their models require a freely 
parametrized reacceleration of electrons, presumably by the motion of turbulent 
magnetic fields (e.g., Seo & Ptuskin, 1994). Strong & Moskalenko (1998) introduce 
a cutoff Zh for the height above the galactic plane above which cosmic rays freely 
escape. They find an upper limit Zh<12 kpc, on the basis of a fit to the 10 Be/ 9 Be 
ratio observed by Ulysses (Connell 1998). This result is "soft": twice the upper 
limit would still be compatible with the ensemble of data (Lukasiak et al. 1994). 
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Moreover, the galactic CR proton distribution extracted from a fit to EGRET 7- 
ray data, actually favours (Strong & Moskalenko 1998) an ad hoc distribution of 
CR sources that is not as well localized in the disk as the conventional supernova- 
remnant sources are (Webber 1997), even if Zh = 20 kpc or more. This point, and the 
necessity to invoke CR reacceleration, indicate that scale heights of the CR electron 
distribution in excess of the 12 kpc "upper limit" may not be out of the question. 
Our results are optimized by a scale height of roughly 20 kpc. Such a large scale 
height is not in contradiction with radio synchrotron-emission from our galaxy if the 
galactic disk and its magnetic field are embedded in a larger magnetic halo with a 
much weaker field. 

In studying the possibility that the diffuse GBR is not extragalactic, one has 
two choices. The first is to extend to high galactic latitudes the elaborate models 
(with many parameters, reacceleration, and ad hoc modifications of the CR-proton 
and CR-electron energy and source distributions) that have been developed to de- 
scribe the intricate nature of the observations at low galactic latitudes (Strong & 
Moskalenko 1998; Moskalenko & Strong 2000). The second is to adopt our very 
naive set of hypotheses and employ a simple cosmic-ray model with, by conven- 
tional standards, a large scale height for CR-electrons. Models of this type (Dar 
& Plaga 1999), wherein cosmic ray sources are directly injected at high galactic 
latitudes, have actually been proposed^. 

In Section 7 we discuss the magnitude and angular-dependence of the two dom- 
inant contributions to the GBR within our model: inverse Compton scattering of 
galactic CR-electrons off the cosmic background radiation and starlight. In Section 
8 we compute the small additive effect of sunlight, and in Section 9 we estimate the 
contribution from external galaxies, which is also sub-dominant. In Section 10 we 
compare our predictions with the data on the intensity and the angular dependence 
of the GBR. The results are very satisfactory and, within our model, lead to the 
conclusion that the GBR can be dominated by the emission from our own galaxy. 
We summarize our conclusions and predictions in Section 11. 



2 The GBR data 

We call "the GBR" the diffuse emission observed by EGRET by masking the galactic 
plane at latitudes |b| < 10°, as well as the galactic centre at |b| < 30° for longitudes 

1 The injector agents would be highly relativistic jets from the birth of compact objects in 
supernova explosions, leading to a CR population permeating a magnetized region of galactic-halo 
proportions and constituting a putative solution to the problem of the origin of the highest-energy 
cosmic rays, a qualitative description of the nuclear CR spectrum, and a possible explanation of 
jetted gamma-ray bursts. 
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|1| < 40°, and by extrapolating to zero column density, to eliminate the 7r° and 
bremsstrahlung contributions to the observed radiation and to tame the model- 
dependence of the results. Outside the mask, the GBR flux integrated over all 
directions in the observed energy range of 30 MeV to 120 GeV, shown in Fig. 1, is 
well described by a power law: 



dF 7 
dE 



(2.74 ±0.11) x 10" 



E 



MeV 



-2.10±0.03 



cm 



s- 1 si' 1 MeV' 1 



The overall magnitude in Eq. ([!]) is sensitive to the model used to subtract the 
foreground (Sreekumar et al. 1998; Strong et al. 1998), but the spectral index is 
not. The EGRET data are given in Sreekumar et al. (1998) for 36 (b, 1) domains, 
9 values for each half-hemisphere. The spectral index is, within errors, extremely 
directionally uniform, as shown in Fig. 3, where we have plotted the EGRET results 
as functions of 8, the observation angle relative to the direction to the galactic centre 
(cos 9 = cos [b] cos [1]). The normalization is less homogeneous, but in directions well 
above the galactic disk and away from the galactic-centre region it has been found 
to be consistent with a normal distribution around the mean value: thus the claim 
of a possible extragalactic origin (Sreekumar et al. 1998). 

In Fig. 4 we have plotted, as a function of 6, the EGRET GBR counting-rate 
above 100 MeV. This figure clearly shows, in three out of the four quarters of the 
celestial sphere, an increase of the counting rate towards the galactic centre. How 
significant is this effect? Let x 2 = X 2 /d-o.f. be the "reduced" x 2 P er degree of 
freedom. The x 2 value for constant flux is 2.6: very unsatisfactory. A best fit of the 
form F = F + Fx (1 — cosO) yields x 2 = 1-3, a very large amelioration (for higher 
polynomials in cos# the higher-order coefficients are compatible with zero: the fit 
does not significantly improve). Note also that at angles with cos# larger than its 
mean value (cos9) =0.0246 {6 < 88.6°), 10 out of the 12 data points are above the 
average flux, while at angles with 9 > 88.6°, 18 out of the 24 data points are below 
the average. The probability for a uniform distribution to produce this large or 
larger a fluctuation is 1.5 x 10~ 4 . 



Even in directions pointing to the galactic disk and the galactic centre, EGRET 
data on 7-rays above 1 GeV show an excess over the expectation from galactic 
cosmic-ray production of 7r°'s (Pohl & Esposito 1998). Electron bremsstrahlung in 
gas is not the source of the 1-30 MeV inner-Galaxy 7-rays observed by COMPTEL 
(Strong et al. 1997), since their galactic latitude distribution is broader than that 
of the gas. These findings also imply that inverse Compton scattering may be much 
more important than previously believed (Strong & Moskalenko 1998; Moskalenko 
and Strong, 2000; Dar et al. 1999). 
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3 The CR data 



The cosmic ray nuclei have a power-law spectral flux dF / dE oc E _/3 with an index 
/3 that changes at two break-point energies. In the interval 10 10 eV < E < Ek ne e 
~ 3 x 10 15 eV, protons constitute ~ 96% of the CRs at fixed energy per nucleon, 
and their flux is (Berezinskii et al. 1990, and references therein): 



dF D 
J - 1-8 


E 


dE 


GeV 


knee < E < E an ki e 



-2.70±0.05 



cm" 2 s- 1 si" 1 GeV" 1 . 



(2) 



In the interval 
to p 2 ~ 3.0, flattening again to fi 3 



3x 10 18 eV, the spectrum steepens from /3i ~ 
- 2.5 above E ankle . 
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The CR flux of electrons (Prince 1979; Nishimura et al. 1980; Tang 1984; Golden 
et al. 1984; Evenson & Meyers 1984; Golden et al. 1994; Ferrando et al. 1996; Barwick 
et al. 1998; Wiebel-Sooth & Biermann 1998), shown in Fig. 2, is well fitted, from 
E ~ 10 GeV to ~ 2 TeV by: 



dFe 
dE 



(2.5 ±0.5) x 10 5 



E 



1 -3.2±0.10 



MeV 



cm 



2 s- 1 sr- 1 MeV -1 . 



(3) 



The terrestrial and solar magnetic fields and the solar wind modify the electron 
spectrum below E ~ 10 GeV, so that the direct observations at those energies may 
deviate from the local interstellar spectral shape. 

Cosmic ray electrons undergo inverse Compton scattering (ICS) off the ambient 
photon baths: starlight and the cosmic background radiation. The spectral indices 
of the GBR and electron spectra can be very simply and successfully related (Dar et 
al. 1999), if the GBR dominantly consists of photons whose energy has been uplifted 
by ICS, as we proceed to show. 



4 The index of the GBR spectrum 

The current temperature, number density and mean energy of the CMB are Tq = 2.728 
K, n « 411 cm" 3 , and e « 2.7 kT « 6.36 x 10" 10 MeV (Mather et al. 1993; Fixsen 
et al. 1996). The galactic starlight (SL) distribution is highly non-uniform, its av- 
erage energy is e* ~ 1 eV. Consider the ICS of high energy electrons on these radia- 
tions. Assume the shape of the electron flux, Eq. (|3|), observed at E > 10 GeV, to 
be representative of the average galactic spectrum. For the energy range of EGRET 
the Thomson limit is accurate even for ICS on SL, and the cross section is 
er T »s 0.65 x 10~ 24 cm 2 . The mean energy E 7 of the upscattered photons, -or AE C , 
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the mean energy loss per collision- is: 



4 / E P x 2 



with €{ = e or 



The ICS photon spectrum originating in our galaxy is the sum of CMB and SL 
contributions: 



dF° 



dF* 



dE dE dE ' K J 

and is a function of the galactic latitude (b) and longitude (1) coordinates. The ICS 
final-photon spectrum -a cumbersome convolution (Felten & Morrison 1966) of a 
CR power spectrum with a photon thermal distribution- can be approximated very 
simply. Using again the index "i" to label the CMB and SL fluxes: 
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where Eg is obtained from Eqs. (J|) by inverting E 7 (e ; ). We postpone to Section 6 
the discussion of the model-dependent normalization factors N*(b, 1) and N (b, 1): 
effective column densities resulting from the convolution of the space distribution of 
CR electrons with those of starlight and of the CMB. Introducing the CR-electron 
flux of Eq. (|), of the form dF e /dE = A [E/MeV] _/3c , into Eqs. (|), we obtain: 



dF 7 _ Nj(b,l) <t t A 
dE ~ 2 



4 ei MeV 
3 m^c 4 
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-, <3c + l 
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cx[E] 



-2.10±0.05 



(7) 



In the energy-range of EGRET, the CMB and SL contributions have the same 
spectral index, as do the small sunlight and external-galaxy contributions discussed 
in Sections 8 and 9. 



The photon spectral index of Eqs. (0), which is related to that of the CR-electrons 
through /5 7 = (/5 e + l)/2, coincides with the measured one, Eq. ([[]). The electron 
spectrum of Eq. (Q) describes the data in the range E e > 5 GeV, so that Eq. (0) 
should be valid above E 7 ~ 100 keV, the typical energy of photons up-scattered 
from the CMB. At E 7 > 50 GeV, at the upper end of the EGRET data, o"t in the 
SL contribution should be replaced by the complete Klein-Nishina cross section, 
implying a steepening of the spectrum. The corresponding effect for the CMB 
contribution is at energy above the EGRET energy range. 

In deriving Eqs. (0), we have assumed that the locally- measured slope of Eq. (0) is 
representative of the index of the spectrum of the electrons suffering ICS to produce 
the GBR, wherever they may be. The spectral index of the diffuse GBR observed 
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by EGRET is independent of direction, as shown in Fig. 3. The statistical test for a 
flat distribution is surprisingly good: x 2 ~ 0.5. This is encouraging support for our 
working hypothesis of an electron spectrum with a universal shape, and of a simple 
and dominant mechanism -ICS- to generate the GBR. 



5 The index of the electron spectrum 

To relate the spectra of CR electrons and protons, we need an estimate of the pro- 
tons' spectrum at their source. A source spectrum dF s /dE with index (3 S ~ 2.2 is ob- 
tained from collisionless shock simulations (Bednarz & Ostrowski 1998) or analytical 
estimates of acceleration by relativistic jets (Dar 1998). The CR spectrum of nuclei 
is modulated by their residence time in the Galaxy, r ga i(E). For a steady source 
of CRs the energy dependence of the observed flux is roughly that of r ga idF s /dE. 
Observations of astrophysical and solar plasmas and of nuclear abundances as func- 
tions of energy (e.g. Swordy et al. 1990) indicate that r gal (E) oc E~°- 5±01 , explaining 
+ 0.5 ~ 2.7, as in Eq. ©. 

Practically all CR acceleration mechanisms invoke an ionized medium that is 
swept by a moving magnetic field, such as would be carried by the rarefied plasma in 
a supernova shell (Bhattacharjee & Sigl 2000) or by a 'plasmoid' of jetted ejecta (Dar 
Sz Plaga 1999). The magnetic field acts as a moving 'mirror' that imparts the same 
distribution in velocity, or Lorentz factor 7 = E/m c 2 , to all charged particles. To the 
extent that particle-specific losses (such as synchrotron radiation) can be neglected 
at the acceleration stage, all source fluxes have the same energy-dependence. For 
electrons below the anticipated 'electron's knee' at E e = (m e /m p ) E knee ~ 2 TeV, we 
expect dFg/dE oc E^ 8 , with /3 S ~ 2.2. Confinement effects preserve this equality for 
ultrarelativistic electrons and protons: their behaviour in a magnetic maze is the 
same. But, unlike for hadrons, the 'cooling' time of electrons -that are significantly 
affected by the ambient radiation and magnetic fields- is shorter than their galactic 
confinement time, r ga i(E), above a relatively low energy. This implies that the CR 
electron spectrum is modulated mainly by the ICS, and not by the confinement 
time. 

Electrons lose energy not only by ICS on starlight and the CMB, but also by syn- 
chrotron radiation on magnetic fields. All of these processes are essentially the same: 
scattering off photons, either real or virtual. The energy loss is governed by the rate 
at which a single electron interacts with the ambient electromagnetic fields, weighted 
by the corresponding average energy density: P = a T c [n^e* + n e + B 2 / (87r)]. Let 
R p (an inverse time) be the production rate of CR electrons, assumed to be constant 
(Berezinskii et al. 1990), and let dn^/dE be their source number-density spectrum. 
The actual density dn e /dE in an interval dE about E is continuously replenished and 
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depleted by electrons whose energy is being degraded by interactions. This leads to 
a steady-state situation in which production and losses are in balance. Using Eq. (|J) 
we obtain: 

1 P — (f 2 — ^ - R M 

3 (m P c 2 ) 2 dE dE p dE ' 1 j 



For a relatively uniform galactic CR electron density, Eq. (|8|) also applies to the 
local electron flux dF e ~ (c/47r)dn e . Substitute the spectrum dn|/dE ~ E _/3s into 
the flux version of Eq. (|§|) to obtain: 

dFe = 3m 2 c 4 R dF e (/3s+1) 

dE 4 (ft- 1)P EdE ' v ; 

For electrons with E < (m e /m p ) Ek nee we deduced that (3 S ~ 2.2 . Thus, (3 S + 1 = 3.2, 
in agreement with the data: Eq. @ and Fig. 2. Above the 'electron's knee' at E e ~ 2 
TeV the spectrum should steepen up by A/5 ~ 0.25, like that of CR hadrons (Dar 
1998). The available spectral measurements extend only to E e < 1.5 TeV. 

The energy density in the CMB is n e = 0.24 eV cm -3 , coincidentally similar 
to that in starlight at our location: n*e* ~ 0.22 eV cm" 3 . If the local CR and 
magnetic energy densities are in equipartition, B 2 /(87r) ~ 1 eV cm -3 , again in the 
same ballpark. The cooling time of electrons in the ensemble of these fields is: 

T (E) ~ -0.22 x 

The galactic escape time of GeV electrons, which should be similar to that of CR 
protons Tg a i(E) oc E~°- 5±01 (Swordy et al. 1990), has a weaker energy dependence 
than that of r cool . At sufficiently low energy, then, r gal < r ool , and processes other 
than Compton- or synchrotron cooling (such as Coulomb scattering, ionization losses 
and bremsstrahlung) become relevant. The slope of Eq. (|9|) should change as the 
energy is lowered. The spectrum of Fig. 2 shows such a change, but it occurs at 
E < 10 GeV, a range in which local modulations would mask the effect. 



E 



GeV 



Gy. 



(10) 



6 The scale height of CR electrons 



The radio emission of galaxies seen edge-on -interpreted as synchrotron radiation 
by electrons on their local magnetic field- offers direct observational evidence for 
CR electrons well above galactic disks (e.g. Duric et al. 1998). For the particularly 
well observed case of NGC 5755, the exponential scale height of the synchrotron 
radiation is 0(4) kpc. If the CRs and the magnetic field energy are in equilibrium, 
they should have similar distributions, and the exponential scale height h e of the 
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electrons ought to be roughly twice that of the synchrotron intensity, which reflects 
the convolution of the electron- and magnetic-field distributions. The inferred value 
h e ~ 8 kpc for NGC 5755 may not be universal for spirals, since h e is very sensitive to 
the density and distribution of CR sources, gas and plasma in each particular galaxy. 
Moreover, the magnetic field may be in equipartition with cosmic rays only where 
the interstellar plasma is dense enough. It is quite possible for the CR electrons to 
be confined in a large magnetic halo with a field much smaller than that in the disk. 
For these reasons we must discuss the observations of our own particular galaxy. 

Traditionally CR electrons and nuclei were assumed to have a distribution that 
snugly fit that of the visible part of the Galaxy -where their conventional sources 
lie- implying a scale height above the plane of the disk of 0(1) kpc (Broadbend et 
al. 1989). As the data and their analysis became more elaborate, scale heights more 
than one order of magnitude larger were discussed (e.g. Strong et al. 1998). Since 
electrons lose energy to the ambient radiation close to their sources, which have 
traditionally been located in the disk, not very well understood CR-reacceleration 
phenomena have had to be invoked (e.g. Seo & Ptuskin 1994). Even with reacceler- 
ation, a conventional distribution of cosmic-ray sources fails to describe the observed 
GBR (Strong k Moskalenko 1998). 

Over the years, Moskalenko, Strong and their collaborators have developed what 
is presumably the most elaborate and detailed understanding of the CR, radio and 
7 observations of our galaxy (Moskalenko et al. 1998; Moskalenko and Strong, 2000; 
Strong and Moskalenko, 1998; Strong et al. 1997; Strong et al. 1998). A crucial 
parameter in their models is the scale Zh of the CR distribution orthogonal to the 
galactic plane, defined as the height above which CRs freely escape, as in a leaky- 
box model. Strong & Moskalenko (1998) conclude that Zh lies between 4 and 12 
kpc. The limits are based on the comparison of the 10 Be/ 9 Be ratio observed by 
Ulysses (Connell 1998) with model predictions as a function of Zh, being all other 
parameters fixed at their adopted values. The dependence of the 10 Be/ 9 Be ratio on 
Zh, shown in Fig. 9 of Strong & Moskalenko (1998) and reproduced here as Fig. 5, 
is very weak for z h > 10 kpc. At z h = 20 kpc, the prediction would be only some 
1.3 standard deviations below the Ulysses central value, and even Zh = 40 would be 
viable: the average of all previous and somewhat less precise observations, compiled 
in Lukasiak et al. (1994) and shown in Fig. 5a, would be in agreement with Zh = 20 
or 40 kpc. For all these reasons and the ones stated in the introduction, we shall 
not refrain from considering scale heights above the 12 kpc upper limit quoted by 
Strong & Moskalenko (1998). 
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7 The CBB and SL contributions to the GBR 



The spectral index of the GBR, derived in Section 4, is independent of the details 
of the spatial distribution of starlight. We have argued that the EGRET GBR data 
support the simple hypothesis of an electron spectral index that is independent of 
location. The predicted GBR index is then also independent of the magnitude of the 
electron spectrum as a function of position. In this section we use a simplified model 
of the electron and starlight distributions to compute the magnitude and angular 
dependence of the CMB and SL contributions to the GBR. 

We adopt h e = 20 kpc (a value obtained from a rough fit of our results to the 
angularly-averaged fluence of the GBR) for the Gaussian scale height of the CR 
electron distribution of our galaxy in the direction perpendicular to the galactic 
plane. For the distribution in p -the radial coordinate orthogonal to the galactic 
axis- we adopt a Gaussian scale height p e = 35 kpc; the results are quite insensitive 
to this parameter. The EGRET GBR data are not precise enough to be "invertible" , 
that is, for the actual high-latitude CR-electron distribution (Gaussian, exponential 
or otherwise) to be disentangled; a fact to be rediscussed anon, in view of our results. 
The distance of the solar system to the galactic centre is d Q ~ 8.5 kpc. The factor 
N o (6',0) in Eq. (Q), which describes the angular dependence of the GBR photons 
due to ICS on the (uniformly distributed) CMB, is: 



N (b,l) = / dr n Exp 
Jo 

h(r, b) = r sinb , 



+ 



do 

Pe 



Exp 



h ( r > b ) V ( P(r,t>,l) x 1 



p(r,b,l) = ([r cos(b) cos(l) -d ] 2 + [r cos (b) sin (l)] 2 ) 1/2 , (11) 
where r is the distance in the direction along the line of sight. 

It is difficult to model in detail the contributionn from ICS on starlight (Hunter 
et al. 1997, Sreekumar et al. 1998). But we are only concerned with this light at 
high galactic latitudes, since the diffuse GBR of interest to us is that measured 
by EGRET by masking the galactic plane and centre. We make a coarse estimate 
by approximating the Galaxy's starlight as that produced by a source at its centre 
with the galactic luminosity L* = 2.3 x 10 10 L Q ~ 5.5 10 55 eV s -1 (Pritchet & van 
den Bergh 1999). The starlight contribution in Eq. (gj) is then of the same form as 
Eq. (|Tl~D, with N traded for N* by the substitution: 

110 ^17^ (r 2 -2rd Q cos(b) cos(l) + d2) " ^ 



For the CMB and starlight contributions to the GBR, averaged over the EGRET 
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unmasked domain, we obtain, by integration of Eqs. @, (|]), (|TT|) , ([12|) : 

-2.10±0.05 



dF 7 
dE 



;2.41±0.55) xlO 



-3 
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MeV 



cnrV^r^MeV- 1 . 
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For scale heights h e and p e similar to the ones adopted (20 and 35 kpc, respectively), 
the CMB and SL contributions are comparable in magnitude, the first scales approx- 
imately linearly with h e while the second is rather insensitive to this parameter. The 
contribution to the CMB from sunlight and external galaxies, discussed in Section 8 
and 9, adds corrections of 6% and ~10% (respectively) to Eq. (|T3"D, the total result 
is shown in Fig. 2. The fitted value of h e is imprecise: the starlight to CMB ratio is 
proportional to e*/eo raised to a very poorly determined power, 0.10 ± 0.05. 

We can use our assumed Gaussian distribution of electrons in a halo, with vertical 
and radial scale heights h e and p e , to compute the diffuse 7-ray luminosity of our 
galaxy, which in our model is dominated by ICS on CMB and SL photons. Using 
Eqs. (H), (|), ([□]), (|3) w e obtain, for the luminosity in 7-rays of energy above E: 

L 7 (>E)^L°(>E)+L*(> E), 
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L°(> E) = 1.31 x 10 



L*(> E) = 3.56 x 10 39 



Pe 



35 kpc 

h e 

20 kpc 



20 kpc 



E 



MeV 



-0.10±0.05 



erg/s, 



1 . 1 + u 

— In 

2 u 1 - u 



E 



MeV 



-0.10±0.05 



erg/s, (14) 



where u = yl — ^1/ pi- A future 7-ray telescope, such as GLAST, could possibly 
see the corresponding glow of Andromeda's halo. 



8 Sunlight contribution to the local GBR 



We are only at a distance 1 Q = 1.5 x 10 13 cm from the sun. This entails a small 
but non-negligible contribution to the locally-observed GBR, resulting from ICS off 
photons in the heliosphere. The corresponding photon flux is described by Eq. ([j]), 
with the substitution of e» by the mean energy e ~ 1.35 eV of solar photons, and 
of Nj by N Q , the solar-photon column density along the line of sight. Let 9 Q be the 
angle between the line of sight and the direction to the sun. Then: 



N (cos# G 



L 



7T 



0- 



sin9 r 



(15) 



For a uniform cos# Q distribution during the EGRET data taking, the average col- 
umn density is N Q = 7rL Q /(16cl e Q ), resulting in a sunlight-induced GBR flux: 

E 



dF° 
dE 



1.32 x 10" 4 



MeV 



cm 2 s 



V^MeV" 1 . 



(16) 
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This contribution is roughly 6% of our galaxy's result, Eq. (]T3|) . At E 7 > 75 GeV, 
the spectrum of Eq. ( |T6| ) should steepen, since ICS should then be described by the 
Klein-Nishina cross section, and not by its low energy Thomson limit. 



9 Extragalactic contribution to the GBR 



To estimate this contribution, some concepts and numbers need to be recalled. 
Hubble's constan' is H = 100 h km s -1 Mpc -1 , with h ~ 0.65; Q m and Q\ are matter 
and vacuum cosmic densities in critical units: O = Q m + y = 1+z is the redshift 
factor. In a Friedman model, the time to redshift relation is dy/dt = —Ho f(y) y, with 
f(y) = [(1 — fl) y 2 + fi m y 3 + f^] 1 / 2 . The luminosity density of the local universe 
(Ellis 1997) is p L = (2.0 ± 0.4) x 10 8 hL Q Mpc- 3 . The combination p L /L* provides 
an estimate of the average number density of 'Milky- Way-equivalent' galaxies. If 
the main sources of CRs are young supernova remnants or gamma-ray bursts, the 
CR production rate ought to be proportional (e.g. Wijers et al. 1997) to the star 
formation rate RsFR,[y], recently measured up to redshift z ~ 4.5 (Steidel et al. 1998). 

The energy of CMB photons up-scattered by electrons at 'epoch y' is proportional 
to T(y) = y T and it is subsequently redshifted by the same factor; hence the 
spectra from distant galaxies should have the same energy dependence as from our 
galaxy. The situation for SL photons is more complicated. Young galaxies are bluer 
than older ones, but this effect is overcompensated by the expansion redshift from 
a relatively low y, onwards. Yet, at the energies observed by EGRET, and for the 
redshift values of 0(1) that dominate the extragalactic contribution, all these blue- 
and red-shifts simply relocate the photon energy, while roughly maintaining the 
slope of the spectrum. For the sum of all galaxies, we estimate: 



dF^ G 1 dL 7 p L c f Rsfr(y) y dy 



dE 4vr EdE L* H J 1 R S fr(0) f(y) y 



r 3 



(17) 



where dL 7 /dE is to be obtained from the luminosity of a Milky- Way-like galaxy, 
Eq. (]I4]). For Rsfr[y] we interpolate the summary values of Steidel et al. (1998). In 
writing Eq. fllTD we have ignored the fact that, above E ~ 10 GeV, absorption by 
e + e~ production on the IR-to-UV background becomes relevant (Salamon & Stecker 
1998), so that the extragalactic contribution should be quenched. 

For Q = Q m = 1 the value of the integral in Eq. ( [I7f) is ~ 0.82; it increases to 
~ 1.08 for a currently more fashionable universe with Q — 1, Q\ = 0.7, Q m = 0.3. 
For the latter case, the result is: 



^p EG r n 1 -2.10±0.05 



7 



dE 



2.48 xlO" 



E 



MeV 



cm'V^r^MeV- 1 , 
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roughly 10% of our galaxy's angularly-averaged result, Eq. dl3|) . 



10 Detailed comparison with the EGRET data 

Our predictions for the magnitude of the GBR and its directional dependence on b 
and 1 are shown in Figs. 5 and 6. In Fig. 5 we display separately the contributions 
from ICS off CMB and SL photons in our galaxy, as well as the uniformly distributed 
sunlight and extragalactic components. In Fig. 7 we compare the total GBR flux: 



obtained by summing Eqs. ©, Eq. flU) and Eq. (|13j), with the EGRET data. 
Our result is a satisfactory fit to the observed magnitude and angular trend of the 
GBR (x 2 = 0.98), a vast improvement over the result for a constant (extragalactic) 
ansatz, for which x 2 = 2-6. Although this agreement would be more meaningful, 
had we used a more realistic model of starlight, a more careful treatment may be 
premature, for the EGRET error bars are large enough to accommodate consider- 
able variations in the input modelling. In a previous analysis (Dar et al. 1999), for 
instance, we obtained a similarly good fit with an assumed constant-density, spher- 
ical CR-electron halo of radius 25 kpc, for which the results have the advantage of 
being simple analytical functions. 

We have neglected various putative extragalactic contributions to the GBR. 
Blazars, because of their beamed emission, may not be very relevant. But CR 
electrons injected directly into intergalactic space by active galactic nuclei, radio 
galaxies or gamma ray bursters, may give rise to a contribution of comparable mag- 
nitude and shape to that of the CR electrons in external galaxies. These or other 
potential sources of GBR photons may imply that our parameters h e and p e have 
been overestimated. But this effect cannot be very large, given our success at de- 
scribing the non-trivial angular dependence of the EGRET data. 



11 Conclusions and predictions 

We have presented a simple understanding of the relation between the spectral in- 
dices of cosmic-ray protons, electrons and the GBR. Accepting the possibility that 
the CR-electron distribution in our galaxy may have a scale height larger than con- 
ventionally believed, we have also argued that the bulk of the GBR could originate 
in our own galaxy. Our modelling is extremely simplistic, but quite successful. 



dF. 



dE 




(19) 
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The predictions specific to our scenario are: 

• The GBR should reflect the asymmetry of our off-centre position in the Galaxy. 

• The halo of Andromeda should shine in gamma rays above a few MeV, with 
a luminosity comparable to that in Eq. ([14]). Likewise, very nearby star-burst 
Galaxies, such as M82, and radio galaxies with large CR production rates, 
such as Cygnus A, may be visible in gamma rays. 

• If the CR-proton and electron acceleration mechanisms are the same, the ex- 
istence of a knee in the observed proton spectrum translates into a related 
result for the power index j3 e of the electron spectrum, which should steepen 
above E « 1.6 TeV by Af3 ~ 1/4. 

• The GBR spectrum should not have the sharp cutoff, above E ~ 100 GeV, 
expected (Salamon & Stecker 1998) for cosmological sources. But it should 
nonetheless steepen around 10-100 GeV, because of the anticipated "knee" 
in the electron spectrum and of the energy- dependence of the Klein-Nishina 
cross section. 

These features of our scenario should be testable when the next generation of 
cosmic-ray and 7-ray satellites (AMS-02 and GLAST) are operational, hopefully by 
2005. In spite of their maturity, cosmic-ray physics and 7-ray astrophysics are still 
young, and thriving. 
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Figure 1: Comparison between the spectrum of the GBR, measured by EGRET 
(Sreekumar et al. 1998), and the prediction for ICS of starlight and the CMB by CR 
electrons. The slope is our central prediction, the normalization is the one obtained 
for h e = 20 kpc, p c = 35 kpc. 
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Figure 2: The primary cosmic-ray electron spectrum (Evenson & Meyers 1984; 
Golden et al. 1994; Ferrando et al. 1996) as measured by Prince 1979 [crosses]; 
Nishimura et al. 1980 [squares]; Tang 1984 [circles]; Golden et al. 1984 [triangles]; 
Barwick et al. 1998 [stars] . The slope is the prediction, the magnitude is normalized 
to the data. 
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Figure 3: EGRET data on the GBR spectral index as a function of 9 the angle 
away from the direction of the galactic center. The line is the predicted spectral 
index. The various plots correspond to the individual half-hemispheres, (a) b > 0, 
1 > 0. (b) b > 0, 1 < 0. (c) b < 0, 1 > 0. (d) b < 0, 1 < 0. 
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Figure 4: EGRET data, organized as in Fig. 3, for the dependence on 9 of the GBR 
intensity above 100 MeV. 
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Figure 5: 10 Be/ 9 Be ratio for the diffusive reacceleration models of Strong and 
Moskalenko (1998). (a) As a function of energy for z h = 1, 2, 3, 4, 5, 10, 15 and 20 
kpc. (b) As a function of at 525 MeV/nucleon, the mean interstellar value for the 
Ulysses data, whose la limits are the dashed lines. The data points in (a) are from 
Lukasiak el al. 1994 (square, Voyagers 1,2; open circle, IMP 7/8; triangle, ISEE 3) 
and Connell 1998 (filled circle, Ulysses). 
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Figure 6: Contributions to the GBR flux above 100 MeV as functions of longitude 1, 
at fixed latitude b, from ICS of starlight (dotted), and CMB in our galaxy (dashed); 
from the total ICS from external galaxies (continuous), and from sunlight (dot- 
dashed). The vertical scale is 10 4 times the number of photons/ [cm 2 ssr]. The 
results are for h e = 20 kpc, p c = 35 kpc. 
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Figure 7: The flux of GBR photons above 100 MeV: comparison between EGRET 
data and our model for h e = 20 kpc, p e = 35 kpc, as functions of latitude at various 
fixed longitudes. The grey domain is EGRET's mask. 
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